Quantum Teleportation of Dynamics and Effective Interactions Between Remote 
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Most protocols for Quantum Information Processing consist of a series of quantum gates, which 
are applied sequentially. In contrast, interactions, for example between matter and fields, as well 
as measurements such as homodyne detection of light, are typically continuous in time. We show 
how the ability to perform quantum operations continuously and deterministically can be leveraged 
for inducing non-local dynamics between two separate parties. We introduce a scheme for the 
engineering of an interaction between two remote systems and present a protocol which induces 
a dynamics in one of the parties, which is controlled by the other one. Both schemes apply to 
continuous variable systems, run continuously in time and are based on real-time feedback. 



Most protocols in Quantum Information Science (QIS) 
are discrete in the sense that they consist of a sequence 
of unitary operations and measurements. Schemes for 
quantum teleportation or dense coding are typical exam- 
ples. These elementary protocols are the building blocks 
of other applications such as quantum repeaters or quan- 
tum computing. However, some implementations are in- 
trinsically continuous. The most prominent example are 
atomic ensembles interacting with light, where schemes 
based on the continuous detection of quadrature opera- 
tors are realized [TH3]. In this system, protocols can be 
performed that are intrinsically deternrinistic and con- 
tinuous in time. Here we address the question how this 
property can be exploited by designing primitives that 
take advantage of this fact. Continuous schemes have 
been devised in several subfields of QIS, for example for 
phase estimation [H |S] , error correction [S] and for the 
preparation and protection of quantum states [7HIU] in 
particular in the context of dissipative schemes [TTMl7| . 
In this letter, we introduce two protocols which achieve a 
qualitatively new goal - to control and transmit quantum 
evolutions between remote locations. We present two re- 
lated protocols, where we consider two remote systems 
which evolve according to a given local dynamics. The 
systems can not interact directly, but they can exchange 
quantum states and classical information. One scheme 
implements an effective non-local dynamics, where the 
two remote systems evolve as if they were interacting 
with each other. The other protocol realizes the quan- 
tum teleportation of a time evolution, which uses the dy- 
namics of one system to steer the evolution of the other. 

The setup under consideration is shown in Fig. [Tk. 
Two spin ensembles interact with a propagating light 
field, which is constantly measured. By performing real- 
time feedback on both samples, one can establish an 
effective interaction between the two separate systems. 
Since this is done continuously, the dynamics of the two 
systems corresponds to the evolution under the desired 
interaction Hamiltonian at any instant of time. Remark- 
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FIG. 1: Dynamical teleportation and creation of an effec- 
tive interaction between two remote systems, a) The Setup 
consists of two atomic ensembles in constant magnetic fields 
oriented along x. Light propagating along z interacts with 
both samples and is continuously measured. The result is 
fed back instantaneously to the atoms, b) Illustration of the 
interaction between the atomic system and the light field in 
terms of discretized spatially localized light modes. 



ably, this scheme results ideally in a joint unitary evolu- 
tion of the two remote atomic systems. This is the case 
even though the protocol is based on measurements yield- 
ing random outcomes and therefore random projections 
of the states involved (see also [T5j). We show that us- 
ing a quantum nondemolition (QND) interaction between 
spins and light, any Hamiltonian, which is quadratic in 
the atomic operators can be realized by tuning the feed- 
back operation only, i.e. the scheme does not require the 
variation of the system parameters. In the ideal case, any 
quadratic Hamiltonian can be implemented perfectly. 

The second protocol realizes a continuous quantum 
teleportation. Quantum teleportation [19] offers a prac- 
tical solution to the delicate task of transmitting quan- 
tum states [Sl| . It is a prerequisite for quantum net- 
works [211 122] and a primitive for quantum computa- 



tion [53]. A standard teleportation scheme consists of 
three separate steps, which involve (i) establishing a 
highly entangled link between the sender and the re- 
ceiver, (ii) a projective measurement, which destroys the 
state to be teleported, and (iii) the recovery of the state 
on the receiver's side by applying a feedback operation. 
Here, we consider a continuous process involving only 
weakly entangled states and measurements which dis- 
turb the quantum state only little at each instant of 
time. In contrast to previous approaches, which transmit 
a (static) quantum state by means of a single feedback 
operation at the end of the protocol, we consider the 
transmission of a whole time evolution using real-time 
feedback. To illustrate this point, we consider the contin- 
uous teleportation between two parties, Bob and Charlie, 
when a time dependent magnetic field is applied to Char- 
lie's system. The resulting displacements on Charlie's 
side translate into a corresponding evolution on Bob's 
side, whose system evolves as if it was placed in a time 
dependent magnetic field. More specifically, it evolves 
as if interacting with a field whose time dependence is 
determined by Charlie's evolution. 

We consider two separate systems that are character- 
ized in terms of continuous variables xj, p\ and %, pu, 
which commute canonically [x,p] = i. We assume that 
both systems can be rotated locally and interact with a 
propagating auxiliary bosonic system via a QND interac- 
tion (Eq. (|2|). For concreteness, we consider two atomic 
spin ensembles interacting with coherent light [21] ■ The 
free evolution of the atomic system is governed by 



H A = 



Oj 



Pi) 



n n 



P 2 u) 



which describes the atomic rotation (Larmor precession) 
in homogeneous magnetic fields along x with Larmor fre- 
quencies fli and fin (see Fig. [1^). In the following, we 
use transformed atomic variables 



\Pi/uJ 
#(Oi/ii) 



= fl(fii,n) 



( Xi/ii 



coaiCljijit) — sin(n I /i I 4) 
sin(Qi/nt) cos(flimt) 



Xl/U 

Pi/ii 



(1) 



which rotate at the Larmor frequencies of the respective 
fields. The light propagates along z, passing both ensem- 
bles. We adopt here a discretized one-dimensional model. 
The interaction time T is divided into N infinitesimally 
small time steps of length r and the incoming light is ac- 
cordingly described in terms of N short pulse pieces. The 
quadratures associated with the n th localized light mode 
are denoted by XL,n an d PL,n with [xL,n>,PL,n'] = i^,,,,,-- 
The n th infinitesimal pulse piece interacts with the atoms 
during the time window from t = (n— l)r to t = w (see 
Fig. [lb) according to the interaction Hamiltonian [25] 



H, 



QND,n 



Nt 



(pi(nr) + piiinr)) Ph,n, (2) 



where k is a dimensionless coupling constant. Such a 
QND interaction can for example be realized using a 



Faraday interaction in atomic vapors [23 I2Z] or in op- 
tomechanical systems [28]. The total Hamiltonian is 

given by H = H-l+H A+Z)n=i ^qnd,„, where H h accounts 
for the free propagation of the light. By virtue of the con- 
secutive QND-interactions shown in Fig. [Ik, the quantum 
state of the atoms is mapped to the x-quadrature of the 
light field. The corresponding discrete input-output re- 
lation for an infinitesimal time step is given by 



-out 
x L,n 



V^(nr)R A ([n~l}r), 



(3) 



with 



Vl^) = (-sia(fl I t),cos(Qit)-sm(Q n t),cos(Cl II t)) , 
Rl(t) = (xi(t),pi(t),xu(t),pii(t)), 

where Or <C 1 has been assumed [29 . The p-quadrature 
is conserved p£ u * = PlV The atomic evolution during an 
infinitesimal interval r is given by 



■-i/iiV 



Pi/n(nT) 



x 1/u ([n-l]r) 

Pi/ii([Vi-1]t) 



K /cos(fi I/n nr) N 



We assume that xl is continuously measured and that 
the result is instantaneously fed back onto the atoms by 
performing a conditional displacement operation. We ap- 
ply here feedback operations with temporally modulated 
gain factors -j^g x j/n(t) and -j^g p j/n(t) such that 



xf u {nr) 



x l/u (nT) 
Pi/u(nr) 



1 (g x j/ji(nr)\ r t (4) 



We show now how this setup can be used to establish 
an arbitrary quadratic interaction between two ensem- 
bles. Using suitable local rotations [30], any interaction 
Hamiltonian for two continuous variable systems, which 
is quadratic in the system operators, can be expressed as 



H = 7 (fJ,H A + isHp) = 7 Zxip 2 



1 



7Pi%2 



(5) 



where \x = \ [Z + ^j and v — \ [Z — -|) [3T]. 7 char- 
acterizes the overall coupling strength of the interaction 
while Z parametrizes the imbalance between the active 
(entanglement creating) and the passive (energy conserv- 
ing) components, Ha and Hp [32] ■ After the light-matter 
interactions, the measured quadrature Xl contains infor- 
mation on both ensembles (see Eq. (pi above). Feedback 
of xl according to Eq. Q , leads therefore to terms which 
correspond to the evolution under both, local and inter- 
action Hamiltonians. The feedback of information of each 
ensemble onto itself leads to an effective evolution accord- 
ing to local squeezing Hamiltonians. In order to suppress 
these contributions, we use the fact that the quadratures 
of each ensemble are mapped to the light field with an 
oscillatory time dependence. For implementing a purely 
non-local evolution, the gain function for ensemble I (II) 



is chosen to oscillate with fi/j (^i), such that informa- 
tion is transferred with high efficiency between the en- 
sembles, while contributions due to the feedback from 
the samples to themselves are out of phase and aver- 
age out. This becomes apparent by considering the time 
varying gain functions g x j/n{t) = -7=g* /h sin(O a/b t) and 



9p,i/ii 



w 



/JV 
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where the coupling matrix -Mr 1 *^ 2 (t) is defined by 

M ai,a 2 (f\— f~9i sin(Oit)sin(n 2 i) .91 sin(f2it)cos(0 2 £) 
Jw si,S2 W I - 52 cos(fti*)sin(fi 2 t) g 2 cos(fti£)cos(ft 2 £) 



AC.i/ii and A/" p .i/ii are noise terms due to the mapping of 
the input light field onto the atomic systems, 



K, l/n (t] 

KanM 



^_/ / 3a/ b sin(n Vb t) «;cos(O vri t)A fx h (ct,0) 
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FIG. 2: Fidelity F versus squeezing parameter r of the light 
field for R — 3 and SI2 = 2fii . The three curves corre- 
spond in ascending order to Q\At — 100 (blue dash-dotted 
line), ili At = 300 (dashed green line) and il\At = 500 (solid 
red line). The gain is chosen such that the resulting interac- 
tion strength corresponds to the underlying QND interaction 
(g — k — 1). The main panel and the inset show the attain- 
able fidelity for realizing a two mode squeezing and a beam 
splitter interaction respectively. Any time evolution under 
a quadratic Hamiltonian can be realized by combining these 
two interactions with fast local rotations [S2, S3 . 



We use here continuous light modes with quadratures 
x L (ct,0) = x^J^ and p L (ct,0) = p^JVr M- As 
shown in the Supplemental Material (SM), the differ- 
ential equations above can be approximated by their 
coarse-grained version using coarse-graining time inter- 
vals At 3> Qjhj, |Oi — fij-j| -1 . In this limit, the cou- 
pling matrices M^ 1 * lead to a negligible contribution 
for Qi y^ J7 2 , since their matrix elements average out. 
Similarly, they can be approximated by a constant diag- 
onal matrix for fli = 51 2 [53]. The noise terms AC,i/n(£) 
and A/" p .i/n(i) give rise to four modes, which are approx- 
imately independent for At ^> ^7/jj> |^j — f2j/| _1 , and 
can therefore be squeezed simultaneously such that their 
contribution becomes negligible. A detailed analysis is 
provided in [34]. For establishing an effective interaction 



according to Eq. (|5j) 

^a/b = ^II/I, 5a/b = 



with 7 = gf, we consider the case 
—gZ Tl , which leads to 



K9_ ( fJTu 

2T 
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By tuning the feedback parameters g a /b, it is therefore 
possible to realize any time evolution, which corresponds 
to a quadratic interaction Hamiltonian. 
As explained above, the setup presented here can not 
only be applied for facilitating an interaction between 
remote systems: using a modified configuration, a contin- 
uous quantum teleportation can be realized. A teleporta- 
tion scheme [TH] involves three parties, - Alice, Bob and 
Charlie. It allows Alice to teleport an unknown quantum 
state provided by Charlie to Bob. In our case, Charlie's 
state is stored in ensemble II and teleported to ensem- 
ble I, which represents Bob, while the light field plays 



the role of Alice. Step (i) in the standard protocol out- 
lined in the introduction corresponds to the interaction 
between the light field and ensemble I which results in 
an entangled state. The distribution of entanglement be- 
tween the two remote sites is realized by the free propa- 
gation of the photonic field. Step (ii) corresponds to the 
interaction of the light with the ensemble II and the mea- 
surement of the x-quadrature of the transmitted light. 
Step (iii) is implemented in the form of a feedback oper- 
ation realizing a conditional displacement on ensemble I, 
which can be done using magnetic fields. Using this pro- 
tocol, the deterministic teleportation of a quantum state 
between two atomic ensembles has been demonstrated 
recently [35j . We present now an extension to a time- 
continuous operation including real-time feedback, which 
facilitates the teleportation of quantum dynamics. For 
this purpose, we consider the special case of the scheme 
introduced above where f^ = fi n = fl and Feedback is ap- 
plied only to the first ensemble g x j(t) = ~7?f$x sin(fii), 

9p,i(t) = 7^5p cos (^)- With this choice, Eq. (W| yields 



^ n ( t )\_1 M n,n () (x I (t) 
^(t))~T M ~ s ^ [t) {pi(t) 



T *»ffp w I p n (t) 



+ 



1 



g x s'm(£lt) KCOs(flt)\ (x\,[ck, 0) 



j 1 \5 P cos(fii) Ksm(flt)Jypi J (ct, 0) 



(6) 



The second term in Eq. ([6]) allows us to control Bob's 
dynamics using Charlie's system. In order to demon- 
strate the transmission of a time evolution, we add an 
extra Hamiltonian corresponding to a time dependent 
transverse magnetic field on Charlie's side -ff cx tra,n = 
a x (t)xn + a p (t)pu. a x (t), where a p (t) are real time de- 



pendent functions. Accordingly, 



fx u (t) 



dt' 



a p (t') 
-a x (t') 

, fcos(nt') 



Pa. 



(7) 



dt 



PL(ct',0), 



where a x (t) and a p (t) describe the time dependence of 
-ffcxtraji i n the rotating frame (see Eq. (fill). The result- 
ing time evolution on Bob's side can be evaluated by 
inserting Eq. (fij) into Eq. ^. As in the general case 
discussed above, the resulting differential equation can 
be approximated by the coarse-grained equation with 
Mg'-g -> -fl for ^At > f, where At is the coarse 
graining time interval. This way, one obtains 



xf n (t) 
Pi n {t)J 2T, 



t 

dt' 



i) 



-ffx« P (i) 



K 

~2T 



gM n (t)+Pu(t)} 



K 

f 



dt' 



Ph(ct,0) 



-g x cos(Qt)\ 

\ g P sm(nt) ) 

-g x sin(Qt) KCOs(Slt)\ fx\ J {ct 1 0) 
5 P cos(Ot) KS\n(Vtt)J \pl(c£,0) 



(8) 



The first term on the right side is equivalent to an effec- 
tive field experienced by Bob. More specifically, Bob's 
ensemble evolves as if it was placed in a transverse mag- 
netic field, whose time dependence is given by Charlie's 
evolution (see Eq. (fij)). This way, the effect of the mag- 
netic field applied at Charlie's side is teleported to Bob. 
This is possible, since entanglement generation, measure- 
ment and feedback are performed continually. In contrast 
to a traditional teleportation, Charlie's quantum state is 
therefore not destroyed in a single step and subsequently 
restored on Bob's side, but rather transmitted continu- 
ously, which offers the possibility to include the effect of 
a time evolution. 

In the remainder of the article, we discuss imperfec- 
tions. To this end, we consider the case, where real-time 
feedback is applied to both ensembles. The following dis- 
cussion can also be applied to evaluating the added noise 
in the teleportation protocol. For analyzing the perfor- 
mance of the scheme, we introduce a figure of merit which 
quantifies the deviation of the realized time evolution 
from the desired one using the Jamiolkowsi isomorphism 
between quantum maps and states. Both, the imperfect 
real map and the ideal one are transformed into their 
corresponding states and then compared as explained in 
the SM. More specifically, the time evolution Et acting 
on the atomic system is described in terms of an entan- 
gled state of twice the system size that can be used to 
teleport a given input state p; n through St, such that the 
output £t(Aii) is obtained. This entangled state consists 
of two copies of a two mode squeezed state with squeez- 
ing parameter R. For R — > oo, which corresponds to a 
quantum state with infinite energy, any input state can 
be teleportet through £y. For finite R, this holds for a 
restricted set of input states. We start by considering the 
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FIG. 3: Fidelity for optimal squeezing r opt versus e R and 
QiAt for Q 2 = 2fii. As in Fig. [5} g — n = l. Fhis plot shows 
the fidelity for realizing a two mode squeezing interaction. 
The results for a beamsplitter interaction are provided in |34j . 



fidelity for fixed R. As explained in the SM, the scheme 
involves two types of imperfections, fast rotating terms in 
the evolution and light noise added to the atomic system. 
If a time window At 3> ^1,^2, |^i — f^l is considered, 
the former are negligible and the latter can be suppressed 
using squeezed input light fields. Fig. [2] shows the at- 
tainable fidelity for fixed R, fli and SI2 = 2f2i versus the 
squeezing of the light field r for different values of Ox At. 
For increasing CliAt, the optimal squeezing parameter 
r op t increases, which leads to an increased accuracy of 
the scheme. For fiiAt — > 00, r opt — > 00 and F — > I. For 
increasing values of the parameter R, correspondingly 
high values of At are required to obtain a good fidelity. 
The required value of At increases with e R as shown in 
Fig. [3] which displays the fidelity F for fixed Larmor fre- 
quencies f^i and Jl 2 = 2^!. This graph features a fine 
substructure since for time windows At = 27r/Oi, lo- 
cal maxima are obtained, which gives rise to stripe-like 
regions of high fidelity. In general, the realization of a 
desired time evolution with a good temporal resolution 
requires high Larmor frequencies. For fixed J7i and O2, 
the attainable precision depends on the temporal resolu- 
tion as shown in Fig.[2]and Fig. [3] We remark that higher 
fidelities can be obtained if a stroboscopic interaction is 
implemented where points of interest in time are chosen 
to coincide with the local maxima. 

In conclusion, we have shown how a time evolution 
can be teleported by means of a dynamical teleporta- 
tion scheme, where entangling operations, measurement 
and feedback are performed continuously and simulta- 
neously. Moreover, we demonstrated how a generalized 
version of this protocol can be used to implement arbi- 
trary quadratic interactions between two separate sys- 
tems. In particular, important tools for QIS, such as a 
purely active (entangling) or a purely passive (quantum 
state transfer) interaction can be obtained. 
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Supplemental Material 

In the following, we discuss imperfections of the proposed scheme for inducing non-local dynamics and show that 
the desired time evolution can be realized perfectly in the ideal case. In Sec. [T] we consider the different types of 
imperfection and explain in which parameter regime the desired time evolution can be faithfully implemented. In 
Sec. [2] we show how the deviation of the implemented time evolution from the desired one can be quantified and 
introduce an adequate figure of merit. In Sec. [31 we include squeezed light in the consideration and calculate the 
resulting fidelity of the atomic evolution numerically. We consider here the general scheme involving feedback on 
both ensembles for generating an effective interaction between two atomic ensembles. The noise contribution in the 
teleportation scheme can be analyzed along the same lines [Si] as the added noise discussed in the following (in 
particular in Sec.[3[). 

1. Imperfections of the scheme 

Throughout the Supplemental Material, we consider two different time scales 

r < fi" 1 < At, 

and assume that the frequencies f^i ,ri 2 , and |f2i — figl take values of the same order of magnitude, which defines a 
time scale f2 _1 . r <C Q^ 1 is an infinitesimally short time interval which will be used to express a continuous time 
evolution in a discretized form. At 3> fl^ 1 is length of the considered time window. 

We consider the setup shown in Fig. la in the main text for establishing an effective interaction between two 
atomic ensembles. The time evolution of the atomic system during an infinitesimal time step of duration r is given 
by 




(1 + G(nr)r) 



Xi(nr) 
pi(nr) 
X\i{ut) 
Pii (nr) 
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N pu n 

N, 



(S.l) 



V2.1l 



The second term on the right is an undesired contribution due to the mapping of the input-quadratures of the light 
field onto the atomic ensembles, 



N x 
N, 



P2,n 



N 



cos(Oit) 
sin(fii£) 

COs(fl2t) 

sin(f22*) 



PL, 



9_ 

N 



( \ sin(fi 2 i) \ 
Zcos(fi 2 t) 
Zsin(fiit) 

y \ cos(O x i) J 



(S.2) 



The matrix G is given by 

/ ^ sm(Q,2t) sin(Jlit) 



G(t) = 



CJH 
Y 



Zcos(£l2t) sin(Oii) 
Z sin(Oif) sin(J7it) 



- ^ sin(fi2t) cos(fiit) \ sin(J72i) sin(J72t) 
-Z cos(ri2i) cos(fiii) Z cos(il2t) sm(£l2t) 
-Zsm(£lit) cos(Oit) Z sin(fiii) sin(il 2 t) 



-^ sin(fl2t) cos(fl2t) \ 
-Z cos(£l2t) cos(5l2t) 
-Z sm(flit) cos(il2t) 



\ j? cos(r2it) sin(£7i£) — j? cos(fl2t) cos(£lit) ^ cos(il2t) sin(£l2^) — j? cos(fl2t) cos(H,2t) J 



The dynamics described by Eq. (S.l) has to be compared to the desired time evolution for the atomic system, 



where X — (xi,pi,xii,pu) T and 



X(t) = G X(t) 
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The scheme involves two types of imperfections. Firstly, the proposed scheme does not implement exactly the atomic 
time evolution matrix Go but rather the time dependent matrix G(t). Secondly, light noise is added to the system (see 



Eq. (S.2|). We start by discussing the first type of imperfection. To this end, we are considering the time evolution 
during a time window At in the absence of added noise. The solution of the differential equation X(t) = G(t)X(t) 
can be approximated by X(t) — e Got X(0) in the limit 

ei,i = (r> 1 At)- 1 «l, 
e 2 . 2 = (r> 2 At)- 1 «l, 
ci,2 = (|^ 1 -fi 2 |At)- 1 «l. 

This is due to the fact that the matrix elements involving two different trigonometric functions average out upon 
integration, for example 

t+At -, pt+M 

dt' sin(fi^') cos(fM') = 0(e iti ), -r- / di'sin(fiji') sin(^t') = 0(e itj ), iori^j, 



AtJt ,--.„,—,-,„, _ v ..,.„ A ^ 



while 



— / dt'sin 2 (^t') = _+0(c M ), — / dt'cos 2 (r!/) = -+0(e M ). 



The effective time evolution is valid on a coarse-grained time scale Ai 3> O x x , S^ x 2 , | f^i — Sl 2 | 1 . This can be 
understood in analogy to familiar examples of effective 2 nd order Hamiltonians such as AC stark shift Hamiltonians. 

As a next step, we include the noise terms in the discussion and consider the complete time evolution 
X{t) = G{t)X(t) + N{t) and its solution 

X(t)=Tetf tG <-^ dT X{0) + N iut , (S.3) 

where T is the time ordering operator and 

pAt 

N [nt = / Te^ tG{r ' )dT ' 'N(r')dr', 



(i 



The components of the integrated noise operator 7V lnt = ( Njgf Nff N^ N^ ) are mutually independent 

[A^ nt ,7Vj nt ]=0, 

for e — ¥ 0, where e = max(ei ! i, e 22 , ei i2 )- Therefore, all four noise contributions can be squeezed simultaneously in 
this limit. Note that the two types of imperfections that occur in this scheme are intimately linked. For e — > 0, the 
first part of the differential equation corresponds to a unitary time evolution, which implies that the noise terms can 
be squeezed such that their contribution to the dynamics becomes negligible. 

2. Comparison of time evolutions 

In order to asses the performance of the scheme, we need to quantify how close the established time evolution is 
to the desired one. We do this in terms of covariance matrices using the Jamiolkowski isomorphism as explained in 



Sec. 2 a The corresponding figure of merit is introduced in Sec. 2 b 



a. Characterization of linear time evolutions using The Jamiolkowski isomorphism 

The Jamiolkowski isomorphism allows one to describe a time evolution in terms of states. Gaussian quantum states, 
i.e. states with a Gaussian Wigner function, are completely characterized by their first and second moments. They can 
be conveniently described in terms of their displacement vector D, with components Di = {Xi), and their covariance 
matrix T, with T^- = ({Xi, Xj} + ), where {•, -} + is the anticommutator |S4| . In the following, we use the Jamiolkowski 
isomorphism in order to characterize a time evolution of the atomic system in terms of a quantum state (T,D), as 



s 



a T 



-*■ £r(Pin) 



:} 



measurement 



Supplementary figure S.l: Characterization of a completely positive map et in terms of a quantum state using the Jamiolkowski 
isomorphism. An operation er acting on quantum system can be represented by the density matrix pj, which corresponds to 
twice the system size, pj can be interpreted as two entangled copies of a quantum state, where er has been applied to one 
subsystem. Given the state pj, the map ex can be (probabilistically) implemented by teleporting an input state pm through 
£t- This can be done by performing a suitable joined measurement on p; n and one subsystem. 



explained in |S5j . To this end, we consider two identical two mode squeezed states with squeezing parameter R. The 
corresponding covariance matrix Tj for the two entangled copies of the atomic system is given by 



/ A(R) C(R) \ 
\ C(R) T A(R) ) ' 



where A(R) and C(R) are 4x4 matrices, 



A(R) = cosh(iZ)l, C{R) = sinh(i?)diag(l, -1, 1, -1). 

In the limit R — > oo, the two mode squeezed state is the improper maximally entangled state |\&tmss,oo) oc X^^Lo l n ' n )- 
A time evolution of the atomic system X(t) = S(t)X(0) + n(t) (where S(i) is a time evolution matrix and n(i) added 
noise) can be characterized by applying the corresponding completely positive map £t to one part of the entangled 
composite quantum state introduced above, which results in the covariance matrix 



J[S,n] 



( SA(R)S T + r noiso SC(R) 



V {DC{R)Y 



A{R) 



nn T + n T n. Fig. S.l illustrates an intuitive explanation of this characterization. The time evolution 



where T DOis 

£t can be applied to an input state p- m by teleporting this state by means of a joint measurement with one part of 
the composite entangled system (see [S5] for a more detailed explanation). For fi->oo, any input state p[ n can be 
probabilistically teleportet in this way, resulting in the output £t{P\u)- For finite values of R, this reasoning can be 
applied to a restricted set of input states. 



A noisefree time evolution et translates into a Jamiolkowski covariance matrix Tj that corresponds to a pure 
state. Deviations of the implemented time evolution from the desired unitary one translate directly into noise added 
to the Jamiolkowski matrix. 



b. Figure of merit 



In order to assess the deviation of the covariance matrix that corresponds to the time evolution implemented by 
the proposed scheme T& tl from the ideal (desired) one, T^, we introduce the figure of merit 



E = maxg 



CHY At -Too) C 

ctr At c 



which quantifies the error. Maximization over all complex vectors C yields 

E = l-m(T E ), 
where 771.(1^) is the minimum eigenvalue of the matrix 



(S.4) 



1 



1 



(S.5) 



This can be seen by noting that 

_, 1 . I &T^C ' 

hi = 1 — mm^; -^ =* 

C [C^T At C_ 

and considering the derivative of the expression A^/Aa* with respect to C, where 

Aa* = C 'TAtC, Aoo = C'TooC. 

The resulting condition for extremality 

AooTAtC ~~ AAtTooC = 
translates directly into 

— -d = T E d, 

Aoo 



where Te is given by Eq. (S.5I and d = v / TaTC 



3. Fidelity of the scheme for squeezed input states 

In this section, we show that the desired time evolution can be implemented perfectly if squeezed light fields are 
used. To this end, we explain how squeezed photonic input states can be included in the derivation and calculate the 
attainable fidelities. 



The time evolution of the atomic system is given by Eq. (S.3|. The total (integrated) noise contribution 
N int = (JV£f , N p f , N£$ , N p f) T can be written in the form 

N int = TY, 

where Y — (xl,1j^l,2) ■■■iXl,n,P'l,i,P'l,2, ■■■,Pl,n) T is the input light field vector and T is a 4 x 2N matrix. We can 
define canonical partners N' = (N' 1 ^ , N n pl , N'™ 2 , N n p2 ) T and normalize the operators such that 

[iVTSiV'f] = iSij, i,j = 1,2,3,4. 

Using these definitions, we can define the vector AT| nt = (n^, N p f , N}$, N$, N'™£, N' p n *, N'^ 2 \ N'%f\ , where the 
subscript indicates that this vector includes the canonical partners. The vector iVg nt can be written in matrix form 

A^ nt = T 8 Y. 

The noise operators contribution to the atomic evolution Njff, NS*, N^, iVS* are independent for e — ► 0. In this 
limit, the total noise commutation relation 



is given by 



[N^,N^]=iS iij 



4 u 
-1 4 4 



where O4 and I4 are the 4x4 zero and unit matrix respectively. For finite values of e, 

&At =T 8 a 00 T s = (Too +0(e). 

Due to the small correction on the right side, the noise modes can not be squeezed simultaneously. We consider 
here a squeezing operation with squeezing parameter r on the input light field which is applied prior to the actual 
protocol. More specifically, we consider the operation on the input light field that corresponds to a squeezing of the 
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Supplementary figure S.2: Fidelity for optimal squeezing parameter r opt of the input light field versus e R and fiiAi for fixed 
Larmor frequencies f2i and f?2 = 2fii. As in Fig. 2 in the main text, the gain is chosen such that the resulting interaction 
strength corresponds to the underlying QND interaction (g — k — 1). This plots shows the attainable fidelity for realizing a 
beamsplitter interaction. 



integrated noise modes |S6j . The fact that a At ^ Coo leads to a tradeoff that prevents the fidelity of the scheme to 
grow continuously with increasing squeezing parameter r. Therefore, there is an optimal squeezing parameter r opt 
for a given value of e. This can be understood by expressing the integrated noise modes in terms of canonically 
commuting modes n$ — (ni,n2,n 3 ,n4,n' 1 ,n' 2 ,n 3 ,n' 4 ) T (with [ni,rij] = i(o'oo)ij), 
the atomic covariance matrix can be written as 



such that the noise contribution to 



-1 noise — & J- n ~r P *■ n' s 

with real coefficients a and /?, where /? is on the order of e. r n and T n > are the covariance matrices corresponding 
to n — (ni,n2,ri3,n4) T and n' = (n' 1 ,n 2 ,n 3 ,n' 4: ) T respectively. Considering a squeezed state of the modes n with 
r„ = e~ r l4 and IV = e r l4 leads to the tradeoff described above. Fig. 2 in the main text shows the fidelity 
F = 1 — E (where E is the error as defined in Sec. |2b| ) versus the input squeezing r for fixed values of Qi, tt 2 and 
R. We analyze here the realization of an active (two mode squeezing) and a passive (beamsplitter like) interaction, 
since any evolution under a quadratic Hamiltonian can be realized by combining these two interactions with local 
rotations |S2( IS3] . The plot shows how decreasing values of e lead to an increased optimal squeezing parameter r opt 
and therefore to an increased accuracy of the scheme. 

For high values of the parame ter R , correspondingly high values of At are required to obtain a good fidelity according 
to the defin ition given in Sec 2 b More specifically, the optimal fidelity depends on the product e R AT~~ 1 as shown 
in Fig. S.2 which depicts F = 1 — E versus e R and At for fixed Larmor frequencies fii and SI2 = 2fii. Apart from 
that, the graph features a fine substructure of local maxima for At = 27r/£7 1 . 



In summary, high fidelities can be obtained for e <C 1, which requires At 3> £l\,£l%, |fii — f^. As outlined 
in the main text, this relation can be understood in terms of a time coarse graining with a coarse graining time 
interval At. The realization of a desired time evolution with a very good temporal resolution requires high Larmor 
frequencies. For fixed Larmor frequencies VL\ and fi2> the attainable precision depends on the temporal resolution 
as shown in the figures. We remark that higher fidelities can be obtained if one realizes a stroboscopic interaction 
where points of interest in time are chosen to coincide with the phase matching condition. 



[SI] In the teleportation scheme, light noise is added to the atomic system (which is described by the terms in the second and 
third line in Eq. (8) in the main text). The corresponding integrated noise operators commute in the limit QAt (where 
At is the considered interaction time) if the gain factors gain factors g x — —g p are chosen. In this case, the desired time 
evolution is teleported perfectly with a phase shift, which is not relevant. The discussion in Sec. [3] can be directly applied 
to the added noise in the teleporation scheme. 
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